Fluctuation pressure of a fluid membrane between walls through six loops 
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I. INTRODUCTION 

Membranes are frequent structures in chemical and bi- 
ological systems. Their dynamic behavior at finite tem- 
perature is of great interest, since their dominant repul- 
sive force is given by thermal out-of-plane fluctuations 
[1, 2]. If the temperature is sufficiently high, the details 
of the potential that inhibits their mutual penetration or 
that causes them to be confined to a certain geometrical 
region are unimportant. Then the membranes' thermal 
fluctuations may be described by a two-dimensional field 
theory with a hard-wall potential that describes their 
mutual interactions and the boundary conditions of the 
space accessible to them. 

In an important class of membranes, their constituent 
molecules are able to move freely within them. The 
thermal fluctuations of these "fluid" membranes are con- 
trolled by their bending rigidity k. The curvature energy 
of such a membrane is, in the harmonic approximation, 
described by 



limit A — ► co. By scaling analysis, the fluctuation pres- 
sure of the membrane has the form [1] 



E=^J A <Px[d\{x)]\ 



(1) 



where the subscript refers to a plane with an area A that 
serves to parametrize the membranes' surface, and where 
ip(x) describes the location of the membrane orthogonal 
to the point x on this plane. For the harmonic approx- 
imation to be valid, the membrane must not fluctuate 
too wildly and thus the temperature must also not be 
too high. It is difficult to describe the membranes' fluc- 
tuations outside the range of validity of the harmonic 
approximation, since then, e.g., overhangs with respect 
to any given plane and steric self-interactions of the mem- 
brane are possible. 

There have been various theoretical approaches to 
compute the pressure of a single membrane between walls 
[1, 3-8] or of a stack of membranes [1, 3-5, 9]. Here we 
consider a fluid membrane between two rigid walls and 
ask what pressure its classical statistical bending fluctu- 
ations exert on the walls. The plane parametrizing the 
membrane is taken to be midway between the enclosing 
walls, which are a distance d apart, and we consider the 
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and we are interested in the numerical value of a. Esti- 
mates of a range from crude theoretical estimates a « 
0.0242 by Helfrich [1] and a w 0.0625 by Janke and Klein- 
ert [3] (this reference also contains an early Monte Carlo 
result a = 0.060 ± 0.003) through Monte Carlo results 



a = 0.079 ±0.002 



by Janke, Kleinert, and Meinhart [4] and 



a = 0.0798 ± 0.0003 



(3) 



(4) 



by Gompper and Kroll [5], and a theoretical estimate 
a w 0.0771 by Kleinert [6] based on the analogy with 
a quantum mechanical particle in a box to a theoretical 
estimate a w 0.0797 by Bachmann, Kleinert, and Pelster 
[7] using variational perturbation theory. Recently, we 
have extended the four-loop calculation in [7] to five loops 
[8] and found a value a ~ 0.0820, outside the error bars 
of the Monte Carlo results. We were, however, unable 
to quote an error bar for our own result. In this work, 
we extend our computation through six loops. Together 
with improved resummation methods, this allows us to 
confirm the disagreement with the Monte Carlo results 
and put stringent error bars around our result. 

Our work is structured as follows. In Sec. II, we briefly 
remind the reader how the hard walls may be modeled 
using an analytic potential and how a perturbative series 
for a may be derived. In Sec. Ill, the central results of the 
technically similar quantum mechanics (QM) problem of 
a particle in a box are listed since they are instrumen- 
tal for extracting a for the membrane problem from its 
perturbative expansion in Sec. V. In Sec. IV, we directly 
resum the perturbative series for a. In Sec. V, we adjust 
the potential modeling the boundary conditions for the 
membrane problem so that the perturbative series for a 
of the QM problem is obtained. Appropriate resumma- 
tion schemes let us then infer the distance of the walls 
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described by the resulting potential, and this information 
is trivially translated into a value of a for the membrane 
problem. In Sec. VI, we summarize and briefly discuss 
our results. 



II. MODELING OF THE HARD WALLS 

Consider a tcnsionlcss membrane between two large 
flat parallel walls of area A separated by a distance d, 
whose curvature energy is given by (1). The d-dependent 
part fd of the free energy density of the system at tem- 
perature T is given by the functional integral 



exp 



Afd\ 
k B Tj 



n 

X 



-d/2 



d/2 



dp(x) exp 



E 



The pressure is then obtained as 



P= — 



djd 
dd 



(5) 



(6) 



and has the form (2) [1, 3]. Our goal is to find the nu- 
merical value of the constant a. 

Following an idea introduced in [6] and utilized also in 
[7, 8], we implement the restriction —d/2 < ip < d/2 by 
adding a potential term m 4 d 2 J d 2 xV(<p/d) to E, where 
V is an even function that is analytic inside a circle 
with radius 1/2 and has sufficiently strong singularities 
at ±1/2. We then expand the potential V in a Taylor 
series in ip and drop the restriction on ip. At the end of 
the calculation, we let m — > to recover the hard-wall 
limit. 

Since the functional form of p in terms of k, d, and T 
is known and since we differentiate fd only with respect 
to d, we set fc^T = n = 1 in the sequel. The energy 
functional may then be written as 



E = 



2 p(x)} 2 + -m^xf 



+ m 4 e d 2 + m 4 e 2k d 2 ^p(x) 2k , (7) 

k=2 ' 

where the e 2 k are the expansion coefficients of the poten- 
tial V. 

The above procedure defines a finite-m version /<j(ra) 
of the free energy density fd of (5), such that fd = 
lim m _>o fd{m)- fd{m) may be expanded in a perturba- 
tive series in terms of vacuum diagrams — i.e., Feynman 
diagrams without external legs [6-8]. The technical de- 
tails of this procedure are described in [8] , and deviations 
from the treatment in [8] are delegated to the appendix. 
The result is that an expansion of fd(m) through L loops 
has the form 



fd(m) 



1 L 



d 2 



(8) 



with the expansion parameter 



1 



m 



■d 2 ' 



(9) 



The perturbative coefficients a\ are functions of the e 2 k- 
Combining (2), (6), and (8), we obtain a finite-g version 
a(g) of a such that an expansion of a(g) through L loops 
has the form 



1 

y 1=0 

of which we need to extract the limit 



a = lim a(g). 



(10) 



(11) 



In Sees. IV and V, we consider several resummation 
schemes for extracting the value of a from a limited num- 
ber of coefficients a;. 



III. QM PARTICLE IN A BOX 

A one-dimensional problem similar to the two-dimen- 
sional case above is finding the ground state energy of a 
QM particle in a one-dimensional box [6, 10] (which, in 
turn, is equivalent to finding the classical partition func- 
tion of a string with tension between one-dimensional 
walls [6, 10, 11]). Introduction of a potential to model 
the hard walls leads to a quantity a(g) parametrizing the 
ground state energy of a particle moving in this potential 
(see [6, 8, 10] and the appendix for details; our notation 
follows [8]). This quantity has a loop expansion of the 
form (10) and due to the trivial topologies of the Feyn- 
man diagrams through two loops, the coefficients a , a\, 
and <22 are identical to those of the membrane case. 

For the particular potential 



Vc(z) = 



1 



2-7T 2 cos 2 (7rz) 



the exact ground state energy is known 
and translates into 



see. 



(12) 



e.g., [10]) 



128 



16 1 4 
+ 2 + ~g 



ttV 
64 



giving the limit 



a = — = 0.07710628438 . 
128 



(13) 



(14) 



for g — > oo. The QM result (13) will be utilized by the 
resummation schemes of Sec. V to extract a for the mem- 
brane problem. 

In the sequel, we will always contrast the membrane 
results with those for the QM problem for the same re- 
summation scheme. 
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TABLE I: Expansion coefficients for the poten- 
tial V c and for a potential V m b that gives the 
QM coefficients a; also for the membrane problem. 



TABLE III: a from VPT as applied in Sec. IV for both the 
QM and the membrane problem. a m b for L = 2,3,4 and 
L = 5 were already obtained in [7] and [8], respectively. 





V c 


v mh 


L 


Qqm 


ttmb 


eo 


1/(2*0 = 


0.0506606 


same 


2 


0.0385531 


0.0385531 


£2 


1/2 = 

7T 2 /3 = 


0.5 


same 


3 


0.0719411 


0.0737974 


£4 


3.28987 


same 


4 


0.0758821 


0.0794726 


£6 


17tt 4 /90 = 


18.3995 


18.0284 


5 


0.0767518 


0.0813538 


C8 


31tt 6 /315 = 


94.6129 


89.5702 


6 


0.0769910 


0.0820175 


do 


691tt 8 /14175 = 


462.545 


419.568 








ei2 


10922tt 10 /467775 = 


2186.57 


1890.91 









TABLE II: Perturbative expansion coefficients for both 
the QM and the membrane problem for the potential V c - 



QM 





1/2^ = 


0.0506606 


same 


ai 


1/8 = 
tt 2 /64 = 


0.1250000 


same 


0,2 


0.1542126 


same 


a.3 


tt 4 /1024 = 


0.0951261 


0.105998 


0,4 







0.026569 


as 


-7r 8 /262144 = 


-0.0361959 


-0.034229 


as 







-0.083246(13) 



membrane 




IV. a FROM DIRECT RESUMMATION OF a(g) 

Knowing only a few low-order coefficients ai, we are 
looking for the g — > oo limit of the series (10). This 
limit corresponds physically to removing the regulator 
that suppresses fluctuations in the infrared. In the con- 
text of critical phenomena, such series have been success- 
fully resummed using Klcincrt's variational perturbation 
theory (VPT; see [12-14] and Chaps. 5 and 19 of the 
textbooks [15] and [16], respectively; improving pertur- 
bation theory by a variational principle goes back at least 
to [17]). Accurate critical exponents [13, 14, 16] and am- 
plitude ratios [18] have been obtained using VPT. 

In this section, we present the results of applying VPT 
directly to the scries (10) as described in [7] and in Sec. IV 
of [8]. We refer the reader to [8] for the details and just 
mention that we only present the results of the q = 1 
version of VPT, since the results for self-consistent deter- 
mination of q from the series (10) remain too imprecise 
even at the six-loop level. 

For the potential (12), which also plays an important 
role in the resummation variants considered in the sec- 
tions below, the €2k are listed in Table I, and the corre- 
sponding perturbative coefficients ai are listed for both 
the QM and the membrane problem in Table II. The 
membrane's coefficients start deviating from the QM co- 
efficients at the three-loop level. At the beginning, the 
deviation from the particular feature of the QM series 
that even loop orders beyond two loops have zero coeffi- 
cients is small. This gives the membrane's series a struc- 
ture that can be expected to be in a transitional phase 



FIG. 1: Q qm (lower lines) and a mb (upper lines) as 
a function of the number of loops L. The horizon- 
tal line is the exact QM result. The dotted-dashed, 
short-dashed, solid, long-dashed, and dotted lines repre- 
sent data from Tables III, V, VII, VIII, and IX, respec- 
tively. Further explanations are given in the main text. 
At the right, the Monte Carlo results (3) and (4) (boxes) 
as well as our final result (27) (diamond) are displayed. 



towards its high-order behavior. However, for the VPT 
resummation scheme to work well and give trustworthy 
results, it is important that the truncated series to be re- 
summed resembles already the behavior at high orders. 
Consequently, the dependence on the variational parame- 
ter in VPT, when applied directly to the series (10), docs 
not develop increasingly flatter plateaus through the or- 
ders considered. Such plateaus are, however, an internal 
consistency check of the method, and we therefore de- 
velop other resummation variants for obtaining a in the 
sections below. Nevertheless, we provide in Table III the 
extension to six loops of Eq. (24) in [8] for the potential 
V c and for comparison also list the corresponding QM 
results, taken from Table I in [8]. The results are also 
plotted in Fig. 1 (dotted-dashed lines), and in spite of 
the above critical remarks they agree perfectly well with 
the results of the more refined resummation variants to 
be discussed below. 

In [19], an attempt was made to extract a from the a/ 
through six loops using so-called factor and root approx- 
imants. However, the achieved accuracy was not high 
enough for any decision about a discrepancy with the 
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TABLE IV: Expansion coefficients of (16) 
for the quantities l/\fV~ c and l/VVmb- 



QM 



"0 


V2ty = 


4.44288 


same 


«2 




-21.9247 


same 


Va 


tt 5 /12V2 = 
-7r 7 /360\/2 = 


18.0324 


same 


v 6 


-5.93242 


10.3394 


v a 


tt 9 /20160V2 = 


1.04555 


-18.7293 


V10 


-tt 11 /1814400v / 2 = 


-0.114657 


-2.24970 


Vl2 


tt 13 /239500800a/2 = 


0.00857287 


25.4600 



membrane 



TABLE V: Results for a for both the QM and 
the membrane problem using the simple resum- 
mation scheme from the beginning of Sec. V. 



L 


Qlqra 


Omb 


1 


0.0625000000 


same 


2 


0.0792468245 


same 


3 


0.0770188844 


0.0845718 


4 


0.0771087134 


0.0817113 


5 


0.0771062388 


0.0816335 


6 


0.0771062850 


0.0818696(2) 



Monte Carlo results (3) and (4). 



V. a FROM ZERO OF POTENTIAL 

Instead of resumming (10) directly, we apply here the 
strategy of Sec. VI of [8]. That is, as a first step we fix 
the e 2 k in (7) order by order such that the expansion of 
a(g) for the membrane problem V m b is identical to that 
of the QM case with a potential V c . The resulting e 2 k are 
listed in Table I. The second step is then to ask where 
the resulting potential V m b(z) has the singularities ±zq 
closest to the origin on the real axis. The scaling relation 
/ cx 1/d 2 when m 2 — allows us then to recover a for 
the membrane case through 



a mb = 4:Z 2 a q 



(15) 



with a qm from (14). Since the nearest singularities of 
V c are of quadratic type, we may assume that the re- 
sulting membrane potential V m b has approximately such 
a behavior. We may therefore assume that l/\/Vmb has 
approximate linear behavior at its first zero — i.e., at ±zq. 

The simplest investigation of z is to truncate the ex- 
pansion of 1/W at L loops, 



(16) 



1=0 



and subsequently numerically determine the first zero of 
the right-hand side of (16). The vu for the QM and mem- 
brane cases are listed in Table IV. The resulting values of 
a may be found in Table V. While the correct QM value 
(14) is approached exponentially fast, the convergence in 
the membrane case is also remarkable. The results arc 
plotted as short-dashed lines in Fig. 1. 

We interpret the fact that the last two differences 
among the membrane values are comparable as a signal 
that the maximum achievable accuracy with the current 
method has been reached. A more refined approach is 
then needed to take into account a likely more compli- 
cated analytic structure of the potential V^b- Let us 
therefore employ VPT to improve the naive resumma- 
tion above. Consider the quantity y/V{z) — y/V(0). At 



TABLE VI: Expansion coefficients of (18) for 
both the QM and the membrane problem. 



QM 



Wl 


2V2/n = 


0.900316 


same 


U'2 


-10/3 = 


-3.33333 


same 


U3 


128^71-/45 = 
-104tt 2 /21 = 


12.6375 


13.1791 


U.4 


-48.878 


-54.6843 


U5 


6904v / 27r 3 /1575 = 
-817847T 4 / 10395 = 


192.214 


235.065 


U6 


-766.379 


-1037.10 



membrane 



least in QM, the singularities of this quantity nearest to 
the origin are simple poles. The resulting series 



2 = 1 



F(z) ee y/V{z) _ JVW) = E 
may then be inverted to 

oo 

z 2 = ][> F '- 



(17) 



(18) 



l=i 



The first few coefficients ui for both the QM potential V c 
and the resulting membrane potential V m b arc listed in 
Table VI. We are interested in finding 



z 2 = lim z 2 {F). 

F^oo 



(19) 



Motivated by the successes of such an ansatz in crit- 
ical phenomena, we assume that the function F can be 
expanded around its first singularities ±z as 



F = Y,u k {z 2 

fe=0 



2\-q/2+k 



4) 



where q = 2 for QM. Inversion of (20) gives 

oo 

- 2 = E <^- 2ro/9 > 



(20) 



(21) 



m=0 



with z 2 — u' q . We may either set q = 2 as in QM in the 
hope that the deviation for the membrane case is small, 
or determine q self-consistently. We use both approaches 
below. 
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TABLE VII: Results for a when the a ; are fixed to 
be those of the QM problem and q = 2 is assumed. 



L 






2 


0.0750000 


0.0750000 


3 


0.0766754 


0.0791616 


4 


0.0769828 


0.0806435 


5 


0.0770794 


0.0812768 


6 


0.0770973 


0.0815743 



TABLE VIII: Results for q and a when the ai 
are fixed to be those of the QM problem and 
q is determined from its own resummed series. 



L 


9qm 


Q!qm 


<7mb 




3 


2.09487 


0.0786643 


2.26290 


0.0856888 


4 


2.05356 


0.0777648 


2.21951 


0.0846057 


5 


2.02049 


0.0772965 


2.11817 


0.0829441 


6 


2.00822 


0.0771659 


2.08666(1) 


0.0825299(2) 



Now apply VPT [15, 16]. In a truncated expansion 



(22) 



2=1 




-lg/2 



= (tF) 1 {l+t 



2/q 



-lq/2 



(23) 



reexpand the resulting expression in t through t L , set 
t = 1, and then optimize the resulting expression in F, 
where optimizing refers to finding appropriate stationary 
or turning points according to the principle of minimal 
sensitivity [20]. That is, we replace 



L-l 



^E 



k=0 



-lq/2 
k 



2/q 



(24) 



and optimize the resulting expression in F. In the limit 
F — > oo of interest to us, this amounts to 



opt. 



E^E 



U=i 



fc=0 



k ) 



(-1)* 



(25) 



which is the L-loop approximation to z 2 — i.e., using the 
expansion coefficients through ul- It turns out that, 
through the order we are working, there is exactly one 
extremum for even L and exactly one turning point and 
no extremum for odd L. This makes the choice of the 
optimization unique at each order. The value of a is in 
each case obtained through (15). 

The results for q — 2 arc summarized in Table VII and 
plotted as solid lines in Fig. 1. The correct QM value (14) 
is approached exponentially fast. The convergence in the 
membrane case is also remarkable. Though the values are 
slightly lower than those reported in Tables III and V, 
they clearly point towards a value of a above the results 
(3) and (4). 

If we refrain from making assumptions about q for 
z 2 (F), we can determine it sclf-consistently by treating 



TABLE IX: Results for q and a when the 
ai are fixed to be those of the QM problem 
and q is determined from optimized plateaus. 



L 


<?qm 


Q!qm 


9mb 




3 


2.09730 


0.0787132 


2.28225 


0.0861317 


4 


2.04990 


0.0777101 


2.21532 


0.0845332 


5 


2.02405 


0.0773337 


2.13061 


0.0830987 


6 


2.00948 


0.0771767 


2.09303(1) 


0.0825984(2) 



first d\nz 2 /dh\F in VPT [13, 16], since it has the same 
q as z 2 (F) and since 



d\nz 2 

lim — = 

F^oo d\nF 



(26) 



by the assumption of a singularity of the potential. That 
is, we resum the expansion of d In z 2 /d In F as detailed 
above and tune q such that optimization with respect to 
F leads to (26). Through two loops, the expansion of 
d\nz 2 /dlnF is q- independent, and we start with L = 3. 
It turns out that through the order we are working, we 
must use turning points for even L and maxima for odd 
L when determining q. For subsequently determining Zq, 
the situation is reverse — namely, as above for q = 2. 

The results for q and a through six loops are listed in 
Table VIII. The results for a are plotted as long-dashed 
lines in Fig. 1. Note how q approaches 2 rapidly for the 
QM problem and that also for the membrane problem a 
value around 2 appears to be approached. 

An alternative to using (26) for the determination of 
q is to tune q such that the plateaus at which the result 
depends least on variations of F are optimized [21]. This 
strategy has been successfully applied in [21, 22] in the 
context of critical phenomena. In practice, this means 
finding F and q such that first and second derivatives of 
the right-hand side of (25) with respect to F vanish (for 
turning points) or such that first and third derivatives 
with respect to F vanish (for extrema). 

The results for q and a through six loops are listed 
in Table IX and are very similar to those of Table VIII. 
They are plotted as dotted lines in Fig. 1. 
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VI. SUMMARY AND DISCUSSION 

In summary, we have computed the constant a 
parametrizing the pressure law (2) of an infinitely ex- 
tended fluid membrane between two parallel hard walls. 
The hard wall was replaced by a smooth potential, al- 
lowing for a perturbative loop expansion for a. Several 
resummation schemes were used to extract the hard- wall 
limit from expansion coefficients through six loops with 
results listed in Tables III, V, VII, VIII, and IX and 
plotted in Fig. 1. 

The values from Table VII on the one hand and Ta- 
bles VIII and IX on the other hand approach each other 
with increasing numbers of loops from below and above, 
respectively A conservative procedure for combining our 
results for a is to average the lowest and highest six-loop 
values for a obtained above (i.e., the values for a from 
Tables VII and IX, respectively) and take their difference 
to be the full error bar. This provides our final result 

a = 0.0821 ±0.0005, (27) 

displayed in Fig. 1. It lies above the Monte Carlo results 
(3) and (4), also displayed in Fig. 1. The simplest expla- 
nation we have to offer for this discrepancy is that their 
error bars, in particular that of (4), may have been cho- 
sen too optimistic. Also, finite-size or other systematic 
effects may not have been taken into account properly. 

On the other hand, it is possible that the treatment 
in our work is inflicted by systematic errors. Experience 
with VPT tells that the internal consistency checks, es- 
pecially the development of increasingly flatter plateaus 
in the optimization procedure with higher orders, are re- 
liable indicators of VPT to work. These checks have suc- 
cessfully been implemented for the procedures of Sec. V. 
Nevertheless, it cannot be ruled out that the method used 
here is not flexible enough to adequately take into ac- 
count the unknown true analytical structure of a(g). 

Another concern is the value of q m b- While the treat- 
ments leading to the a values listed in Tables VIII and 
IX and thus leading to the upper dotted and long-dashed 
curves in Fig. 1 determine q m b self-consistently, the QM- 
inspired value g m b = 2 was somewhat arbitrarily chosen 
to obtain the values for a in Table VII. What if the value 
of q m b describing a(g) best differs from 2? A slightly 
larger value (but below the ones from the self-consistent 
determinations of q m b) leads to slightly increased values 
of a and therefore also to a larger mean value and smaller 
error bar in (27). On the other hand, a slightly smaller 
value leads to slightly decreased values of a and therefore 
also to a smaller mean value and larger error bar in (27) . 
In both cases, the q m b = 2 curve may still approach the 
correct result. For an optimal value of q m b below 2 this 
may happen by decreasing a values at higher orders (that 
such a behavior is possible in principle is easily tested by 
setting q m b to a value above 2.085, which causes the cor- 
responding six-loop value for a to drop below that at five 
loops). The smooth behavior and slow flattening of the 
q m b = 2 curve lets us believe, though, that such a drop, 



if present, should be very small. The best g m b value can 
thus be expected to lie at most slightly below 2, leading 
to only a small decrease of the values of Table VII and the 
upper solid curve in Fig. 1 and to a slightly lower mean 
value and larger error bar than given in (27). It is re- 
assuring that the results of both the direct resummation 
of a(g) employed in Sec. IV and the naive resummation 
from the beginning of Sec. V lie very close to the mean 
value of (27). 

All things considered, we are rather confident about 
our result (27), but further studies of the system are re- 
quired to settle the question of the correct value of a. 



APPENDIX: FEYNMAN DIAGRAMS 

In [8] , we have described at length how recursion rela- 
tions along the lines of [23, 24] can be used to construct 
the vacuum diagrams needed for the computation of the 
perturbative coefficients ai in (10) through a given loop 
order, and there is no need to repeat the derivation here. 
Given the Feynman diagrams, a vertex with 2k lines rep- 
resents a factor 

-m 4 d 2 ^- k h 2k , (28) 

and the perturbative coefficients a; are obtained from the 
sum of Hoop diagrams as 

ai = - ^2 ci_ n gi_ n Ii_ n , (29) 

n 

where c/_„ is a combinatorial factor, gi_ n is a monomial in 
the 62k, and //_„ is the corresponding momentum space 
integral. The integration measure is J d D k/(2n) D with 
D = 1 for QM and D = 2 for the membrane. The mem- 
brane propagator carrying a momentum k is given by 
l/(k 4 + m 4 ), while the QM propagator carrying a mo- 
mentum k is given by l/(k 2 + to 2 ). 

From the list of diagrams through five loops in [8] it is 
obvious that a major reason for the rapid increase of the 
number of Feynman diagrams with the number of loops 
is a proliferation of the momentum-independent one-loop 
propagator insertions of the form Q . A simple measure 
to reduce the number of diagrams to be considered is a 
one-loop resummation where we absorb the above inser- 
tion into the parameter to in the propagator and at the 
end reexpand the resulting modified perturbative series 
in powers of g as was carried out in [25] and treated on a 
more formal level in [24]. Note that many other resum- 
mations in the form of momentum-independent propa- 
gator and vertex corrections are possible, leading to a 
further reduction of the number of Feynman diagrams. 
However, at the current level of computing vacuum dia- 
grams through six loops, this is unnecessary. In any case, 
the diagrams most difficult to evaluate are always those 
with the full loop topology and remain after any such 
resummation. 
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To implement the one-loop resummation in the mem- 
brane case, we must compute diagrams with a modified 
propagator H such that 

-(4) 



G 12 



H 12 



12^1234^34, 



(30) 

where the notation of [8] has been employed. Writing 



G- 1 = k 4 +n 
translates into 

k 4 + r 

or 



and H 



k 4 + M 4 , this condition 



fc 4 + M 4 - |m 4 d- 2 e 4 M- 



M 2 
m 2 



M 2 
m 2 



3 

r 49 



where we have used (9) and the one-loop result 
d 2 k 1 1 



/ 



Define Z(g) by 



(2tt) 2 k 4 + m 4 8m 2 ' 



M 2 = Z(g)m 2 



Although (32) can be solved analytically for Z(g) 
more useful for our purposes to write 



z(g) = i 



and extract 



Cl 



fc=i 



3 C 
464 



and the recursion relation 



Cfc 



3 fe_1 

1 i= i 



fe-2 fc 



(31) 
(32) 

(33) 

(34) 
it is 

(35) 
(36) 
fc > 1. (37) 



Through six loops in the vacuum diagrams, we need Z(g) 
through g 5 and obtain 



z(g) 



1 + 



fe 4 <7- 



27,2 2 
32 € i9 



4- ^f 3 o 3 
t 16 £45 



8505 4 4 
2048 e 4a 



+ Z2£ e ^+0(«,°). 



(38) 

The same resummation (30) can be implemented for 



the QM case. Writing G~ 



and H 



M 2 , the condition (30) translates now into 
fc 2 + m 2 = k 2 + M 2 - 6m 2 d- 2 e 4 M- 



TABLE X: Numbers of diagrams for low loop orders. 



number of loops I 





1 


2 


3 


-1 


5 


6 


7 


diagrams 


1 


1 


1 


3 


7 


24 


83 


376 


diags. after one-loop resum. 


1 


1 


1 


2 


3 


11 


29 


125 


diags. with Hoop topology 


1 


1 





1 


1 


5 


8 


37 



1. Discard all vacuum diagrams with one or more one- 
loop propagator insertions Q , with the exception 
of the two-loop diagram, which changes sign (here, 
the combinatorics do not work out; loosely speak- 
ing, it is undefined which part of the diagram is the 
insertion) . 

2. Compute all remaining diagrams with the replace- 
ment m — > M in the propagators. 

3. Replace M 2 — > Z(g)m 2 for the membrane case and 
M — ► Z(g)m for the QM case and reexpand the 
perturbative series in powers of g. 

In Table X, we give the original numbers of diagrams at 
some low loop orders, the numbers left after our one-loop 
resummation, and the numbers of diagrams with the full 
respective loop topology. The latter is necessarily the 
same both before and after the one-loop resummation. 
In Table XI, we list all diagrams through six loops left 
after the one-loop resummation. Also given are their 
combinatorial factors c;_„ , their coupling constant factors 
gi- n , and the values /;_„ of the corresponding integrals for 
M = 1 for both the QM and the membrane problem (the 
multiplying power of M can immediately be inferred from 
the number of loops and propagators of a given diagram) . 

The techniques for evaluating the integrals are ex- 
plained in [8]. With the exception of the membrane 
integrals have been evaluated to the precision given either 
in momentum space or in both momentum and configura- 
tion space. For Iq.^, the indicated precision could only be 
obtained in configuration space. Since the slightly lower 
precision of 76-5 introduces the main computational error 
into the determination of a at the six-loop level, we have 
indicated the ensuing numerical error in the other tables 
of this work, where applicable. 



or 



M 



rn 



M 



in 



0. 



(39^ABLE XI: Diagrams l-n (nth /-loop diagram) through six loops, their combinato- 
rial factors ci_ n , coupling constant factors gi. n , and values 7j_ n of the corresponding 
integrals for M — 1. D — 1 and D — 2 correspond to the QM and membrane prob- 
lems, respectively. 

(40) 



where we have used g = 4/md 2 [8] and the one-loop result 

50 dk 1 1 
, 27r fc 2 + m 2 2m 



(41) 



This time we define Z(g) = M/m. Then Z(g) is the same 
as in the membrane case considered above. 

The resulting modifications in the Feynman rules are 
for both the membrane and the QM case: 



l — n 


diagram 


Cl-n 


9l-n 


jEt 1 


r D = 2 
'l-n 


0-1 




1 


— «0 


1 


1 


1-1 


O 


1 

2 


1 


-1 


1 

~4 


2-1 




-3 


— 64 


1 

4 


1 

64 


3-1 




12 


4 


1 

32 


4.04576 x 10~ 4 
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diagram 









7 










15 



288 



360 



105 



2592 



10368 



5760 



12960 



1320 



360 



2700 



5040 



945 



3 
512 



1 

64 



1 

16 



5 

4096 



19 



7 

6144 



7 

3072 



3 

1024 



5 

1024 



1 

192 



1 

128 



1 

64 



1 

128 



1 

32 



-n | diagram | ci-n \ 9l-n \ 



1 

512 



1.63237 x 10~ 5 



5.05719 x 10" 5 



1 

4096 



7.55133 x 10~ 7 



1.04187 x 10" 



6.71540 x 10~ 7 



1.50770 x 10" 



2.04047 x 10" 



3.95093 x 10" 



3.76084 x 10" 



6.32149 x 10- 



65536 



6.32149 x 10" 



32768 



6-1 



6-2 



6-3 



6-4 



6-5 



6-6 



6-7 



6-8 



6-10 



6-11 



6-12 



6-13 



6-14 



6-15 




DOC 



7 



^7 










7 






124416 



248832 



497664 



331776 



414720 



207360 



155520 



622080 



155520 



34560 



51840 



35 



131072 



71 



294912 



367 



1179648 



269 



24576 



25 
65536 



65 



147456 



12288 



19 

24576 



5 

8192 



1 

1024 



5 

8192 



5 

4096 



7 

6144 



3.74650 x 10" 



3.12644 x 10" 



4.52023 x 10" 



2.85447 x 10" 



2.37861(5) x 10" 



6.39380 x 10" 



6.27924 x 10" 



5.50218 x 10" 



1.30234 x 10" 



9.43917 x 10" 



8.39425 x 10" 



1.70039 x 10" 



1.02301 x 10~ 7 



2.46933 x 10~ 7 



1.88463 x 10~ 7 
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